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O ■ Abstract 

&\ I We propose a renormalon-inspired resummation of QCD perturbation theory 

based on approximating the renormalization scheme (RS) invariant effective charge 
beta- function coefficients by the portion containing the highest power of &=g(UJV— 
Qh! 2JVy)> for SU(./V) QCD with Nf quark flavours. This can be accomplished using 

exact large-TVy all-orders results. The resulting resummation is RS-invariant and 
the exact next-to-leading order (NLO) and next-to-NLO (NNLO) coefficients in any 



RS are included. This improves on a previously employed naive resummation of the 
leading-6 piece of the perturbative coefficients which is RS-dependent, making its 
comparison with fixed-order perturbative results ambiguous. The RS-invariant re- 
summation is used to assess the reliability of fixed-order perturbation theory for the 
e + e~ i?-ratio, the analogous r-lepton decay ratio R T , and Deep Inelastic Scatter- 
ing (DIS) sum rules, by comparing it with the exact NNLO results in the effective 
charge RS. For the i?-ratio and R T , where large-order perturbative behaviour is 
dominated by a leading ultra-violet renormalon singularity, the comparison indi- 
cates fixed-order perturbation theory to be very reliable. For DIS sum rules, which 
have a leading infra-red renormalon singularity, the performance is rather poor. In 
this way we estimate that at LEP/SLD energies ideal data on the i?-ratio could 
determine a s (Mz) to three-significant figures, and for the R T we estimate a the- 
oretical uncertainty 5a s {m T ) ~ 0.008 corresponding to 5a s (Mz) — 0.001. This 
encouragingly small uncertainty is much less than has recently been deduced from 
comparison with the ambiguous naive resummation. 



1 Introduction 



There has been a great deal of interest recently in the possibility of identifying and resum- 
ming to all-orders the Feynman diagrams which dominate the large-order asymptotics of 
perturbation theory [1-6]. 

More precisely consider some generic dimensionless QCD observable D(Q), dependent 
on the single dimensionful scale Q, with a perturbation series 

D(Q) = a + dia 2 + d 2 a 3 H h d k a k+1 H , (1) 

where a = a s /n is the renormalization group (RG) improved coupling. The perturbative 
coefficients dk can themselves be written as polynomials of degree k in the number of 
quark flavours, Nf] we shall assume massless quarks. 

4 = d [k] N k + dl~ 1] N k f l + • • ■ + df . (2) 

The leading d k coefficient corresponds to the evaluation in each order of perturbation 
theory of a gauge-invariant set of Feynman diagrams containing chains of k fermion bub- 
bles. Techniques for evaluating such diagrams exactly in all-orders have been developed, 
and rather compact results for QCD vacuum polarization |7|, |), Deep Inelastic Scattering 
(DIS) sum rules ||, and heavy quark decay widths and pole masses |ll|] obtained. 



The resummation of such diagrams provides a gauge-invariant effective charge in QED, 
but in QCD one would need to include chains of gluon bubbles and ghosts, and the 
isolation of a gauge-invariant subset of diagrams providing an analogous QCD effective 
charge is problematic [|H|. One knows on the grounds of gauge invariance, however, 
that part of the result should be proportional to b k , where b is the first beta-function 
coefficient, b—^(HN—2Nf), for SU(iV) QCD. Since for large-iV/ one must obtain the 
QED result one can substitute N f =(^N - 36) in the l N f -expansion' of equation (2) to 
obtain a '6-expansion' 

4 = 4 fc v + dt 1] b k - 1 + • ■ ■ + 4 0) , (3) 

where % —(— l/3) k d k , so that exact knowledge of the leading- Nf d k to all-orders implies 
exact knowledge of the leading-6 d k k ^ to all-orders as well. 

The existence of so-called infra-red (IR) and ultra-violet (UV) renormalon singularities 
in the Borel plane at positions z=2l/b with / a positive or negative integer, respectively, 
means that in large-orders the perturbative coefficients are expected to grow as 4 ~ b k k\. 
As shown in reference fl3| this singularity structure leads to the expectation that the 
'leading-6' term when expanded in powers of Nf should, asymptotically, reproduce the 
sub- leading coefficients. That is, expressing d k b k as 

d[% k = d[ k] N k + d[ k - 1] N k ~ l + ■■■ + d[ k - r] N k - r + • • • + df , (4) 



one can show that, 



S k ~ r] - d[ k - r] 



(5) 

so that for fixed r and large k the sub-leading W/-expansion' coefficients are reproduced. 



i + 0( I) 
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As demonstrated in reference [p]] for both the e + e~ Adler .D-function and DIS sum 
rules this asymptotic dominance of the leading-6 term is already apparent in comparisons 
with the exact next-to-leading order (NLO) and next-to-NLO (NNLO) perturbative co- 
efficients, d\ and g?2- For instance for SU(iV) QCD the first two perturbative coefficients 
for the Adler .D-function are given by [O, [15] 



.6557V + , (6) 

d 2 = .086iV 2 + N f (-1A0N - ^) + (2.ION 2 -. 661-^) . (7) 

These are to be compared with the leading-6 terms 

d!\ ] b = .3456= -.115/V/ + .634/V , (8) 
d (2) b 2 = .7766 2 = .086N 2 f - M8N f N + 2MN 2 . (9) 

The subleading N and NfN, N 2 coefficients approximate well in sign and magnitude 
those in the exact expressions in equations (6) and (7). 

Notice that the level of accuracy with which the sub-leading coefficients d k are 
reproduced is far in excess of that to be anticipated from the asymptotic expectation of 
equation (5). This is a rather weak statement which implies only that d£ ^ should be 
reproduced to 0(1/ k) accuracy for fixed r and large k on expanding d^b k , whereas the 
d^k i r= k) coefficient, which is leading in the 1/7V expansion for a large number of colours, 
is reproduced accurately for k=l and k=2. 

The above observations suggest that there should be some merit in resumming to all- 
orders the 'leading-6' terms, even though many features of the approximation remain to 
be clarified. In a number of recent papers [1-6] such a programme has been carried out. 
We shall refer mainly to the results of reference M. 

Following I] we define 

oo 

D^ = j24 L) ct h+1 (10) 

where d^ = d^b k (d^ = 1). We can also consider the complementary sum over the 
sub-leading b terms 

oo 

D {NL) = j2df L) a k+1 (11) 



k=l 



where d[ NL) = d k - d[ L) . Hence 



D = + . 



In reference [0 the summation in equation (10) was defined using Borel summation, 
with the IR renormalon singularities principal value (P.V.) regulated. Whilst such a 
summation can be performed the results for D^ L ' and D^ NL ^ are separately dependent on 
the chosen renormalization scheme (RS), the sum of the two being formally RS-invariant. 
Changing the RS changes the definition of the renormalization group (RG)-improved 
coupling and hence the 'a' appearing in the summation changes. In the full sum this 
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change is precisely compensated by the RG transformations of the coefficient d k under 
the change in RS. However, by restricting oneself to the 'leading-6' piece of the coefficients 
this exact compensation is destroyed and the resulting sum D^ L ' is fatally RS-dependent 
[p], [l|J. One response to this, which has been adopted in references [2-6], is to artificially 
restrict the RG-transformation of 'a' to that contributed by the first term in the beta- 
function only. With this restriction is then RS-invariant since RS changes are exactly 
compensated for at the 'leading-6' level. In our view, however, one can and should do 
much better than this. For several QCD observables one has exact results for the first 
two perturbative coefficients d\ and d%, usually in the modified minimal subtraction (MS) 



scheme with renormalization scale fi=Q [|14| , [lq , |T7| , [18fl . What is clearly needed is a 
resummation in which the exact NLO and NNLO contributions are included, and an 
approximate resummation of the higher orders performed, in such a way that the full sum 
is explicitly RS-invariant under the full QCD RG transformations. In this way one would 
have a test bed for assessing the reliability of fixed-order perturbation theory in any RS 
by seeing how it differed from the RS-invariant resummed result. As one reduced the 
energy scale Q (e.g. the centre of mass-energy in e + e~ annihilation) one could also assess 
how the reliability of fixed-order perturbation theory deteriorates. 

In this paper we shall show in section 2 how such an improved resummation can be 
carried out, and will use it in section 3 to assess the reliability of fixed-order perturbation 
theory for the e + e~ .R-ratio, the analogous r-decay ratio R T , and DIS sum rules. Section 
4 contains discussion of results and section 5 our conclusions. 



2 RS Dependence and RS-Invariants 

We begin by briefly reviewing the RS dependence of the coupling 'a' and the perturbative 
coefficients df.. We refer the reader to reference |]19| for more details. 



The RG improved coupling 'a' satisfies the beta-function equation 
da 



din- 



— a 2 (l + ca + c 2 a 2 + 



+ c k a k H ) = -(3(a) 



(12) 



where r = Mn^, with [i the renormalization scale and A the dimensional transmutation 
mass parameter of QCD. Here 6 and c are universal with 



-{\\C A - 2N f 
8 b 



11 C A C F 



4° A 



(13) 



where for SU(iV) QCD C A =N and C F =(N 2 -1)/2N. 



As shown by Stevenson |2(| the RS is labelled by the parameters r, C2, C3, . . .; the 
higher beta-function coefficients are not universal and characterise the RS. Integrat- 
ing up equation (12) with a suitable choice of boundary condition |2U| one finds that 
a(r, C2, C3, • • • , Cfc) is the solution of the transcendental equation 



1 , ca f a 1 
— hem h / ax 

a 1 + ca Jo 



1 



+ 



1 



x 2 B(x) x 2 (l + cx) 



(14) 
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where B(x) = (1 + cx + c 2 x 2 + c 3 x 3 + • • ■ + CkX k + ••■). 

For consistency of perturbation theory one finds that di(r), d 2 (r, c 2 ), ds(r, c 2 , C3), ■ • •, 
c^(t, c 2 , C3, • • • , Cfe), • ■ •. The combination 

/0o (Q)= r -d 1 ( T ) = Mn| (15) 

is RS-invariant. 

The explicit functional dependence of the dk on RS is conveniently obtained by consid- 
ering the special RS, the effective charge (EC) scheme pi| , in which d%=d 2 =- ■ -=dk=- • -=0, 



so that D=a, and in this scheme the renormalized coupling is the observable itself. From 
equation (15) this RS will correspond to the choice of parameter r=p (ensuring di=0). 
To determine the remaining parameters, C2=P2, C3=P3, ■ ■ •> Ck=Pk, • ■ - , characterizing the 
EC RS, one proceeds as follows. From equation (12) we have that for two RS's, barred 
and unbarred, 

m = , (is) 

where a, a denote the couplings in the respective schemes RS and RS. If the barred RS 
is chosen to be the EC scheme then 

(3(a) = p(a) = a 2 (l + ca + p 2 a 2 + h Pka k + •••)> (17) 

with a=D. Then equation (16) gives 

P(D) = ^{<D)) , (18) 

where a(D) is the inverted perturbation series. 

By expanding both sides of equation (18) as power series in D and equating coefficients 
one obtains 

P2 = c 2 + d 2 — cdi - d\ 

p 3 = c 3 + 2d 3 - 4cM 2 - 2d lP2 - cd\ + 2d\ (19) 



Since po =r — d\ is RS-invariant we can use di itself, rather than r, to label the RS. 
Rearranging equation (19) we can then obtain 

d 2 (di, c 2 ) = d\ + cdi + (p 2 - c 2 ) 

5 1 

d 3 (d u c 2 ,c 3 ) = dl + -cdf + (3p 2 - 2c 2 )d 1 + -(p 3 - c 3 ) (20) 



The result for d n (di, c 2 , ■ ■ ■ , c n ) is a polynomial of degree n in d\ with coefficients involving 
p n ,Pn-i,- ■ ;c and c 2 ,c 3 ,- • -,c n ; such that d n (0,p 2 ,p 3 ,- ■ -,p„)=0. The p 2 ,p 3 ,- • -,Pn,- • ■, are RS- 
invariants which completely characterise the QCD observable D. They are independent 
of the energy scale Q but do depend on the number of active quark flavours, Nf. Given 
just these numbers the perturbative coefficients in any RS can be obtained from equations 



4 



(20). To construct RS-invariant resummations the strategy will be to approximate the 
RS-invariants p k , and then use equations (20) to obtain the approximate perturbative 
coefficients in any arbitrary RS. In this way invariance under the full RG transformations 
of QCD is guaranteed. 

The pk invariants can be organised as an expansion in b, with 

Pk = + P t l) b k - 1 + ■ ■ • + P f + V 1 • (21) 

The b~ l term arises from the fact that in a 'regular' RS such as minimal subtraction the d k 
are polynomials in b of degree k ||22|| , whereas the corresponding beta- function coefficients 
Cfc are polynomials in b of degree k-1 with additional b~ l terms (c.f. the expression for 
c—ci in equation (13)). The RS-invariant combinations in equation (19) in principle could 
contain arbitrary inverse powers of b, but RG considerations guarantee that only b~ l terms 
remain |j22j| . Thus bpk is a polynomial of degree k+1 in b. 

The effective charge beta-function p(D) (equation (17)) will contain Borel plane sin- 
gularities at the same positions as those in D(a) |23| and hence one should expect a 
weak asymptotic result analogous to equation (5), with the p k ^b k+1 term asymptotically 
reproducing the iV/-expansion coefficients of bpk. For the Adler D-function and DIS sum 
rules the level at which this works is again far in excess of that to be anticipated from the 
asymptotic result. The p k term involves only combinations of the d k k \ with for instance 
p!> 2 ' )=c 4 2 ' ) — {d^) 2 , and so the p k can be obtained to all-orders given the exact leading- TVj 
all-orders results. 

For the Adler D-function (D [|IJ) one has the exact result for SU(iV) QCD (where the 
"light-by-light" contribution D is excluded, see [||) 

bp 2 (D) = -0.0243AT 3 . + (0.553N - 0.00151— )N 2 f 

I 

+ (-3.32iV 2 + 0.344 + 0.0612— )N f (22) 

+ (3.79iV 3 - 1.45iV - 0.337^) . 
This is to be compared with the 'leading-6' piece 

b 3 pf\D) = 6 3 (4 2) - (4 1} ) 2 ) = 0.6566 3 

= -0.0243]^ + OAOlNNf - 2.21N 2 N f + 4.04iV 3 . (23) 

Notice the good agreement of the sub-leading NN 2 , N 2 Nf, and iV 3 coefficients. 

For the DIS sum rules (polarized Bjorken or GLS, K [p]]) one has the exact result 

bp 2 {K) = -0.0221M + (0.513iV + 0.00665-biV 2 

j N J 

+ (-3.29iV 2 + 0.505 + 0.0143^)^ (24) 
+ (3.85iV 3 - 1.73iV - 0.337^) , 

which is to be compared with 

b 3 p ( 2 2 \k) = 0.5976 3 

= -0.0221^5 + 0.365iViV 2 - 2.01iV 2 iV7 + 3.68iV 3 , (25) 
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again the NNj, N 2 Nf, and iV 3 coefficients are well reproduced. 

In both cases in the large-iV limit (Nf—0) the RS-invariant p 2 is approximated to 
better than 10% accuracy. The 20% level agreement of the sub-leading coefficients does 
not, unfortunately, guarantee that the overall RS-invariant is reproduced to the same 
accuracy for all N, Nf since there are large numerical cancellations. For instance for 
N=3 and jV /= 5 one has p 2 (L>)=-2.98 (exact), whereas b 2 p ( 2 2) {D) =9. 64. 



3 RS-invariant Leading-6 Resummations 

Approximating the RS-invariants pk by p^ = b k p^ one can now define the RS-invariant 
resummation 

oo 

D (L*) = £ d j>) a *+i , (26) 

where in a general RS d k L *\di, c 2 , C3, • • • , Ck) is obtained by replacing p^ in equations (20) 
by p[ L \ so that 



j(L*) 



1 

d[ L * ] = d 1 

d i 2 L *\d 1 , c 2 ) = d\ + cdi + (p 2 L) - c 2 ) 



(27) 



df*\d u c 2 , c 3 ) = rf 3 + ^crf? + (3p? } - 2c 2 H + l -{pf ] - c 3 ) 



Notice that, unlike the strict 'leading-6' approximation of equation (10), the NLO coef- 
ficient di is now included exactly. If an exact NNLO calculation exists then the exact 
p 2 can be used and P3,P4,- • ■ approximated by Pz\pi\' ' so that d 2 (in any RS) is in- 
cluded exactly. In any case the all-orders sum in equation (26) is formally RS-invariant, 
and can be compared with the NLO, NNLO, N 3 LO,- • -,N"LO,- • ■ fixed-order perturbative 
approximations to assess the accuracy of the fixed-order results, 



Tl 

£(J>)(n)=£ d ( i *) a *+l. 



(28) 



k=0 



The next task is to define the all-orders resummation in equation (26). If we consider 
the EC RS then a=D and 



r = po — 6 In 



_Q_ 



df s (p = Q), 



where for later convenience we have expressed the RS invariant po in terms of the MS 
scheme NLO coefficient with p=Q, c^ s (/i = Q), which is customarily what is computed, 
and where A^jg is the universal dimensional transmutation parameter of QCD. Equation 
(14) in the EC scheme with x 2 B(x)=p(x), the EC beta-function of equation (17), then 
yields 



1 , cD 
D + TTcD 



b\n-^--df s (p = Q)- [ D dx 
A Ms ' Jo 



+ 



p(x) x 2 (l + cx) 



(29) 
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As discussed in reference the EC beta-function p(x) is of fundamental significance 
since 

and so it can be partially reconstructed from measurements of the energy evolution of the 
observable. Given p(x), D(Q) is specified by the solution of the transcendental equation 
(29). The resummation D^ L *^ of equation (26) will correspond to the solution of equation 
(29) with p(x) replaced by p( L *)(x), where 

oo 

p^(x) EE x\l + CX + P2 X 2 + A L) X k ) ■ (31) 

k=3 

For the observables to which we shall apply the resummation exact NNLO results exist 
and so we have included the exact p 2 , rather than p 2 . . 

We can define p( L *) using the principal value (P.V.) regulated Borel sum results for 
D^ L \a) of equation (10) obtained in reference 

D iL \a) = P.V. / dze- z/a B[D (L) ](z) , (32) 
Jo 

where B[D^ L '](z) denotes the Borel transform which potentially contains poles at z=Z[=j- 
(1=1, 2, 3, • • •) corresponding to infra-red renormalons (IR;), and at z=—z\ corresponding 
to ultra-violet renormalons (UVj). The IR; singularities are intimately connected with 
the operator product expansion (OPE) for the observable in question, and the chosen 
regulation of the IR singularities determines the definition of non-perturbative condensates 

In reference |1J results have been derived for the e + e Adler D-function (D) and the 
polarized Bjorken (or GLS) DIS sum rules (K). For these Euclidean quantities one can 
obtain the regulated Borel sum of equation (32) as sums of exponential integral functions 
Ei(Fzi) and Ei(-Fzi), where F = -. The resulting expressions for D^(F) and K^ L \F) 
split into UV and IR contributions are given in equations (48,49) and equations (52,53) 
respectively in reference [J]. Results are also obtained for two Minkowski quantities, the 
e + e _ i?-ratio (R) and the analogous r-decay ratio (R T ). Expressions for ES l \F) and 
R^ L \F) are given in equations (60,64) and equations (69,70) respectively of reference 
in terms of generalized exponential integral functions. 

Equation (18) at the leading-6 level (j3(a)=a 2 ) yields 



p^(x) = (a^(x)f^M 



(33) 



da 

where a^ L \x) is the inverse function to D^ L \a), i.e. D <<L \a <yL \x))=x , and explicitly 

oo 

P W(*)eex 2 (1 + 5>1 Z V). (34) 

k=2 

p( L >(x) can then be straightforwardly obtained from the corresponding D^ L \F) expres- 
sions of |T] for the various observables. For a given x one numerically solves 

D {L) (F(x)) = x , (35) 
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Figure 1(a): NLO and NNLO fixed-order perturbation theory, the naive RS-dependent 
leading-6 resummation, and RS-invariant leading-6 resummation, for R at Q — 91 GeV 
plotted against 'a'; C2 = has been assumed. 



to obtain F(x), and then from equation (33) 



x) = -^ L \F) 



F=F{x) 



Finally comparing equations (31) and (34) one has 

p( L *\x) =p( L \x)+cx 3 + p ( 2 



{NL) x A 



(36) 



(37) 



where p { 2 NL) = p 2 - p 2 L) . 



This p 



(£*), 



x 



can be inserted in equation (29) and the integral 



performed numerically. Given a value of Aj^g and including the known exact NLO result 
for d^ s (p = Q) one can then solve the transcendental equation (29) for D=D^ L *\ Con- 
versely given D=D ( - L *^=Ddata, from the experimental measurement of the observable, one 
can solve equation (29) for A^g. By varying Q, with A^^ and d^ s (p = Q) evaluated with 
the number of active quark flavours, Nf, changing across quark thresholds, one can study 
the Q-dependence of D^ L *\Q). The resummed result 7J)( L *) can also be compared with 
N n LO fixed-order perturbative results. Since d\ and d 2 are exactly included in any RS 
one has 7J( L *)( 1 )= 7J( 1 ) ) £)( L *)( 2 ) = £)( 2 ) ; where 7J)( L *)(™) denotes the truncations of equation 
(28), and denote the exact NLO and NNLO results. 

In Figure 1(a) we have plotted as the dashed curve the leading-6 resummation R^ L \a) 
versus the coupling 'a' for the e + e _ i?-ratio with Q=91 GeV, the t'Hooft scheme corre- 
sponding to B{x)—1 + cx, and minimal subtraction have been assumed with A^=200 
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Figure 1(b): As for Figure 1(a) but for K at Q — 1.5 GeV. 

MeV. There is a monotonic RS-dependence as discussed in reference Jl|. Noting that r 
is related to 'a' using equation (14) one can use 'a' to label the exact NLO and NNLO 
approximants, R^(a), R^(a,c 2 ). The dotted line shows R^(a), and the dashed-dotted 
line gives R^ 2 \a, 0). We have chosen C2=0 to avoid adding an extra axis to the plot. The 
solid line gives the RS-invariant resummation w- L *^> . We note that the fixed-order results 
agree best with the R^ L *^ resummation in the vicinity of the stationary points with respect 
to variation of the RS. This is to be anticipated since the Principal of Minimum Sensitivity 
(PMS) |20j choice of RS avoids the inclusion of potentially large UV logarithms connected 
with the choice of renormalization scale jn|. A similar statement holds for the NLO and 
NNLO results in the EC scheme ||19|| , R^(EC) and R^(EC), corresponding to solutions 
of equation (29) with p( L *^ in equation (31) truncated. These are numerically very close 
to the PMS approximants. The 'optimized' PMS/EC fixed-order NNLO approximant is 
thus seen to be very close to to the RS-invariant resummed result for the i?-ratio at LEP 
energy, indicating that the approximated effect of including N 3 LO and higher corrections 
is small, and thus suggesting that one can in principle accurately determine A^ given 
ideal data. 

In Figures 1(b) and 1(c) the analogous plots for the DIS sum rule K at Q=1.5 GeV, 
and for the r-decay ratio R T (Q=m T =1.78 GeV) have been given. A|j^=320 MeV has 
been assumed. 

In contrast to Figure 1(a) the differences between the fixed-order results and the RS- 
invariant resummations are clearly much larger. Thus at these lower values of Q the 
significance of N 3 LO and higher effects is apparently much greater, and the reliability 
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Figure 1(c): As for Figure 1(a) but for R\ at Q = 1.78 GeV. 



(3) 

with which Ai^ can be determined correspondingly less. We shall quantify this more 
carefully in just a moment. 

In Figure 2(a) we plot, for the e + e" i?-ratio at Q=91 GeV, the fixed-order pertur- 
bative results R^ L *^ (EC) (equation (28)) for n=2 (NNLO) and higher orders (crosses) 
compared with the RS-invariant resummed result M L *^ (dashed line). A^=200 MeV has 
again been assumed. We could of course have chosen to plot the fixed-order approximants 
in any RS, for instance MS with n=Q, but as discussed in connection with Figure 1(a), 
we expect the 'optimized', EC or PMS, choice of RS to approach the resummed result 
more rapidly. We stress that the fixed-order R^ L *^ n \EC) approximants correspond to 
the solutions of equation (29) with p(x) replaced by the truncation of p^ L *\x) in equation 
(31). Since p 2 is included exactly R^ 2 \EC)=R (2 \EC). 

As can be seen from Figure 2(a) the N 3 LO and higher fixed-order results are indis- 
tinguishable from the resummed result with the chosen vertical scale, and there is only a 
small shift between the NNLO and resummed results. Evidently fixed-order perturbation 
theory in the EC scheme seems to be working very well for the i?-ratio at LEP/SLD 
energies. 

In Figure 2(b) we show a similar plot for R at Q=h GeV. A^=279 MeV has been 
assumed. Clearly the approach to the resummed result is somewhat less rapid. The slight 
oscillation of successive fixed-order approximants above and below the resummed result is 
explained by the dominance of the UVi singularity at z=— | in the Borel plane, which is 
the closest to the origin for the i?-ratio. This singularity is responsible for alternating sign 
factorial growth of the perturbative coefficients. Beyond order n=12 the amplitude of the 
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Figure 2(a): Comparison of fixed-order EC perturbation theory (crosses) with the 
RS-invariant resummation (dashed line) for R at Q = 91 GeV. 
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Figure 2(b): As for Figure 2(a) except at Q = 5 GeV. 
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Figure 2(c): As for Figure 2(a) except at Q — 1.5 GeV. 



oscillations increases dramatically, and the fixed-order approximants diverge increasingly 
from the resummation. This is precisely what one would expect to see on comparing the 
Borel sum of an alternating sign asymptotic series with its truncations. We note that 
a similar oscillating behaviour with eventual wild oscillations setting in would also have 
been apparent in Figure 2(a) had we used a finer vertical scale. ^From the large-order 
behaviour one would not expect the wild oscillations to set in until n > 50. 

Figure 2(c) finally shows the corrsponding plot for R at Q—1.5 GeV, with as 
above. The approach to the resummed result is still slower, and the oscillations have only 
just become established when they increase wildly beyond n=9. Nonetheless even at this 
low energy fixed-order perturbation theory is approximating the resummed results, albeit 
much less well. 

This reasonable performance of fixed-order perturbation theory for R is to be con- 
trasted with the situation for the DIS sum rules K. Figures 3(a), 3(b), and 3(c) are 
plotted at the same values of Q as the corresponding Figures 2 for R, and with the same 
vertical scales to enable direct comparisons. In Figure 3(a) at Q=91 GeV we see a much 
slower approach to the resummed result. The fixed-order EC approximations then track 
the resummed result between sixth and tenth order and then for n=12 there is a dramatic 
breakdown. The Borel plane singularities nearest the origin for K are now [|I|] IRi at 
z=t, and UVi at z—— |. It is the presence of the IRi singularity which leads to fixed- 
sign factorial growth of the perturbative coefficients and a consequent deterioration in 
the performance of fixed-order perturbation theory. The relative deterioration compared 
to the i?-ratio can be seen even more clearly in Figure 3(b) at Q=5 GeV. There is a 
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Figure 3(a): Comparison of fixed-order EC perturbation theory (crosses) with the 
RS-invariant resummation (dashed line) for K at Q = 91 GeV. 




Figure 3(b): As for Figure 3(a) except at Q = 5 GeV. 
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Figure 3(c): As for Figure 3(a) except at Q — 1.5 GeV. 

monotonic increase in successive orders with no tendency to track the resummed result. 
The stair-like pattern, with neighbouring odd and even orders roughly similar in low or- 
ders, follows from a partial cancellation between the fixed-sign (IRi) and alternating-sign 
(UVi) behaviours. 

In Figure 3(c) we see that at Q—l. 5 GeV fixed-order EC perturbation theory is a poor 
approximation to the resummed result for the DIS sum rules. Only n=2 and n=3 are 
shown since for n >4 fixed-order perturbation theory is not defined in the EC scheme. If 
p(x) has a zero at x=x* (where x* > is the closest zero to the origin) then equation (29) 
has a solution D = D*, with D* < x*. This will be the case if the expansion coefficients 
of p(x) have alternating factorial behaviour, at least in either odd or even orders. If, 
however, the coefficients have fixed-sign factorial growth, as is the case for the DIS sum 
rules, then p(x) will have no positive zeros. In this case equation (29) may fail to have a 
solution, the condition for this being that in the limit as D — > +oo the right-hand side of 
equation (29) is negative. 

In Figure 4 we show the analogous plot for the hadronic r-decay ratio, R T . Here 
evidently Q — m T — 1.78 GeV, and the same as above has been assumed. Fixed-order 
EC perturbation theory is seen to be working reasonably well with oscillating behaviour 
around the resummed result which becomes wild for n > 5. Notice, however, that the 
performance is much worse than that of R at the comparable Q — 1.5 GeV in Figure 2(c). 

We can summarize the behaviours exhibited in the foregoing figures by plotting the 
energy dependence of R and K. Figure 5(a) shows R^ (solid line), R^(EC) (dotted 
line), and R^(EC) (dashed-dotted line), plotted versus In Q/ GeV over a range equivalent 
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Figure 4: Comparison of fixed-order EC perturbation theory (crosses) with the 
RS-invariant resummation (dashed line) for R r . 



to Q = 1 — 91 GeV. Flavour thresholds in Q at m& = 4.5 GeV, m c = 1.25 GeV, have 
been assumed and values of A^ chosen as above. Figure 5(b) gives a similar plot for 
K. We note the reasonably satisfactory behaviour for R, in particular at all energies the 
NNLO EC approximation is closer to the resummed result than the NLO, as one would 
hope. In contrast for K below Q ~ 3 GeV the NLO becomes closer than NNLO to the 
resummed result, making the use of fixed-order perturbation theory dubious. Notice in 
addition that the vertical scale in Figure 5(b) is much coarser than that in Figure 5(a). 

We would finally like to use the RS-invariant all-orders resummation to assess the likely 
accuracy to which A^jj^ can be determined for various observables at various energies. 

In Table 1 we have given the A^' } values obtained by fitting the NLO, NNLO EC 
fixed-order perturbative and the RS-invariant resummed results to the central values of 
the data for R(Q=91 GeV) |25|, R{Q=9 GeV) |2§, R T , and K{Q 2 =5 GeV 2 ). The value 
for R T is that obtained in B by averaging those obtained from the leptonic branching 



ratio [E7H and r-lifetime measurements 112811 : we have corrected for the small estimated 



power corrections ||. 

For K(Q 2 = 5GeV 2 ) we have taken the GLS sum rule result of the CCFR collaboration 
2£j corrected by subtracting off the Q~ 2 higher-twist corrections suggested by reference 
3(J, so that 

K(Q) = Kccf R (Q) + (Q - 27 ^ 2 °- 14) • (38) 
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Figure 5(a): NLO, NNLO fixed-order results in the EC scheme, and RS-invariant 
mation, for R plotted versus lnQ/GeV over the range 1 < Q < 91 GeV. 
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Figure 5(b): As for Figure 5(a) except for K. 
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Observable 


Energy Q/GeV 


N f 


Data 


A^^/MeV fitted to experiment 










NLO 


NNLO 


Resummed 


R 


91 


5 


0.040 ± 0.004 


9 c 9 +190 


9QO+228 


9Q«+232 




9 


5 


0.073 ± 0.024 


399lfl 


5I6I423 


OO 1 _44 3 


R T 


1.78 


3 


0.205 ±0.006 


3191? 


3871^ 


404±i| 


K 


2.24 


3 


0.154 ±0.075 


437l^ 9 7 


3791*1 


49fi +386 
4ZO -303 



Table 1: Values of adjusted to fit the predictions of NLO, NNLO fixed-order results 
in the EC scheme, and the RS-invariant resummation to the experimental data for R, R T 
and K. 



The errors have been combined in quadrature. 

The results in Table 1 are encouraging in that they indicate rather small differences 
between the NNLO EC and resummed fits. From these differences one would estimate that 
one could determine A^l to an accuracy of ~ ±3 MeV at LEP/SLD energies given ideal 
data for R, corresponding to determining a s (Mz) to three significant figures. Needless to 
say even with ideal data undetermined finite quark mass effects would in fact introduce 
far larger uncertainties. 

At Q = 9 GeV A^l would apparently be determined to an accuracy of ~ ±20 MeV. 
The data for R(Q = 91 GeV) imply NNLO a s (M z ) (MS) values of a s (M z )=0.122±0.012; 
NNLO EC and resummed are the same to the quoted number of significant figures. 

For R T a comparison of the NNLO EC and resummed fits would suggest that A^± 
could be determined with a precision of ~ ±15 MeV. One finds a s (m T )=0.320 ± 0.005 
(NNLO EC) and a s (m r )=0.328 ±0.005 (resummed). Evolving up from N f = 3 to JV> = 5 
assuming; the flavour thresholds noted above yields A^=253lg MeV (NNLO EC) and 

A^l=267±18 MeV (resummed), corresponding to a s (M z )=0.119± 0.001 (NNLO EC) and 
a s (Mz)=0.120±0.001 (resummed). A conservative estimate of the theoretical uncertainty 
is then 5a s {M z )=0.001. 

If taken seriously the above estimate of the accuracy with which A^l (a s (Mz)) can 
be determined from R T measurements is very reassuring, and clearly indicates that this is 
indeed potentially the most reliable determination. The uncertainty is somewhat smaller 
than has been assumed based on more naive estimates of the size of the 0(a 4 ) perturbative 
coefficient [31]. It is much smaller than 6a s (m T ) = 0.05 inferred by Neubert in reference 
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H based on comparison of the exact 0(a^) NNLO perturbative result in the MS scheme 
with fi=m T , and a straightforward resummation of the leading-6 terms, which is essentially 
our (c.f. equation (10)), with a= Qa (™- T ) _ As can be seen from Figure 1(c) the dashed 
curve R^ L \a) lies above the RS- invariant resummation (solid line) for a= " s (™- T ) ~ 0.12, 
and there is a strong 'a' dependence in this region. The NNLO EC result, and indeed 
the MS /i = m T NNLO result, are both much closer to the RS-invariant resummation. 
The implication then is that the rather large difference between the exact fixed-order and 
naive resummed leading-6 results found in || is just a reflection of the inadequacy of the 
naive resummation, which was our original motivation for improving it. 

We finally turn to the GLS sum rule results in Table 1. Whilst the A^ values for NLO, 
NNLO, and resummed are in reasonable agreement with that obtained for R T , we note 
once again the worrying feature that the NLO result is closer to the resummed than the 
NNLO. We have also had to assume and correct for sizeable power corrections, based on 
the modelled suggestion of reference [5(1. Also note the very large errors which reflect the 
difficulty in reconstructing the sum rule by combining data from various DIS experiments 
EjJ. Clearly K will not be competitive with R T as a way of determining Aj^f. 



4 Discussion 

Before giving a summary of our main conclusions we would like to discuss several ways 
in which we could improve or extend the RS-invariant resummations, and mention some 
technical issues related to them. 

The first concerns the analytical continuation between the Euclidean Adler D-function 
and the Minkowski quantities R and R T . This will imply definite relations between the 



corresponding RS-invariant s pj?, pf, pf T . For instance for R one has |32 



P R 2 


= p?- 


1 

12 


pf 


= pf- 


5 

12 


pf 


= pf- 


1 

12 



1.2 2 
71 



9tt 2 (39) 
360 



„D , -7 2\lJ2„2 , 1 l4„4 



The procedure we have used to construct R^*^ involves resumming the effective charge 
beta-function using the exact pf and the leading-6 approximations to pf , k > 2. This 
means the — -^cb 2 n 2 analytical continuation term in pf , or the — j^c 2 b 2 7i 2 in pf, have 
been omitted since they are sub-leading in b. Since p® is known exactly we could also 
improve the resummation by using the exact p® in evaluating the — |p^6 2 7r 2 term in pf. 

One could envisage an improved resummation p^**'* (x) incorporating these extra terms. 

p ( £** ) (x)=p { £*\x)+p R (x), (40) 
where the extra terms to 0(x 7 ) are explicitly 

p R ( X ) = ~^cbW - (|p?< NL) + y 2 C 2 )bW + 0(X 7 ) + ■■■ . (41) 
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This resummation to all-orders can be accomplished in principle by representing R as 
a contour integral involving D | 32fl . Using Z)( L *) in the integrand would formally produce 
ft(L**) correS p 0nc li n g to the above effective charge beta- function p^**^ , but one would have 
to evaluate £H L *) at complex values of Q. Similar remarks apply to R T . One might note 
that in the NNLO case, where we can compare with the exact result, p2^ is only a good 
approximation to the exact p2 for Nf ~ or for large Nf. Hence one could conclude that 
the uncertainties in the basic approximation are such that the attempted improvement is 
not warranted. Nonetheless it would be very worrying if any of our conclusions for R, R T 
changed on including these extra terms. We hope to check this in a future work. 

Another aspect of the resummations which requires elucidation is the way the re- 
summed p( L *' effective charge beta-function is obtained by numerically inverting the P.V. 
regulated Borel integral representation of D^ L \a), as detailed in equations (35)-(37). From 
equation (30) we see that p(D(Q)) is directly related to the Q-evolution of the observable 
D(Q), and is therefore of central physical importance in studying power corrections. One 
might then imagine defining 

p(D) = -Reg." / dz e^ D B[p](z) + p*%(D) , (42) 
Jo 

where B[p] denotes the perturbatively defined Borel transform of p. This will contain 
singularities at the same positions in the Borel plane as D(a) P5| , and to control the 
IRz infra-red renormalon singularities on the positive- z axis the integral will have to be 
regulated, denoted "Reg." . There will be an additional p^p{D) incorporating the power 
corrections (e _1//D terms) whose precise definition will depend on the chosen method of 
regulation p4| . 

If B[p] is defined in the leading-6 approximation we can then ask if the first term in 
equation (42) with P.V. regulation exactly reproduces the p^ L *^ defined by numerically 
inverting the P.V. regulated D^ L \a). 

This can be reduced to a simpler problem. Consider 

D(x) = P.V. r dze- z/x B[D](z) 
Jo 

(43) 

roo 

a(x) = P.V. / dze- z/x B[a](z) 
Jo 

where B[D] denotes the perturbatively defined Borel transform of D(a), and B[a] denotes 
the perturbatively defined Borel transform of the inverse function of D(a), which can 
unambiguously be defined by formally transforming the coefficients of the power series 
D(a). With these definitions one can then ask whether D(a(x)) = x exactly or whether 
there are additional e _1//x terms. Existing results on such problems are in short supply 



34]], but unfinished work in progress |35[ strongly suggests that the relation D(a(x)) = x 
does hold exactly. Unfortunately the result probably only holds for P.V. regulation. The 
pragmatic reason for studying this question is that B[D^ L '](z) is given by rather simple 
expressions as a sum of poles whereas B[p](z) will have an extremely complicated 
form. Hence it is impractical to construct B[p](z) directly, and the numerical inversion 
route is the only possibility. 
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The properties of the function p(x) fix the infra-red properties of the observable. 
For instance if p(D*) = then D — > D* as Q — > 0. It has been suggested [0, [37| 
that such infra-red freezing is supported by a wide body of indirect phenomenological 
evidence. In reference the assumption of universal infra-red behaviour of an effective 
coupling a e ff(k) has been used to interrelate power corrections for different observables. 
It is interesting that for all the observables we have studied in this paper D^ L \a) has a 
maximum value, D max , say. This means that p^ L \x) is undefined for x > D max . If p is to 
be defined in the infra-red this is presumably a signal that power corrections have to be 
included beyond a certain point. An interesting consistency check on this interpretation 
is that if only ultra-violet renormalon singularities are present then D^(a) does not have 
a maximum. In particular if, as is the case for R and R T , the UV singularities are single 
poles, then D^ L \a) increases monotonically and p^ L \x) will be defined for all x. The 
absence of IR renormalons is consistent with there being no power corrections, or at least 
they are not constrained by the large-order perturbative behaviour. We hope to take up 
this question of IR behaviour and constraining the form of power corrections in a later 
work. 

A final underlying issue which needs further clarification is the explanation for the 
excellent performance of the leading-6 approximation itself. For all the cases where exact 
NNLO QCD calculations exist the leading-6 approximation not only gives exact results 
for perturbative coefficients and pk RS-invariants in the large- Nf limit, but remarkably 
is also an excellent (~ 5% level) approximation in the large- N limit of a large number of 
colours. As pointed out unfortunately it may be a rather poor approximation in-between 
these extremes, for Nf=5, N=3, for instance. Although the sub-leading iV/-expansion 
coefficients are reproduced remarkably well (~ 20% level). 

A possible Feynman-diagrammatic explanation runs as follows. In the large- N limit 
of QCD only planar diagrams contribute, 't Hooft has shown that if one restricts oneself 
to UV-finite planar diagrams perturbation theory converges ||39|| . As far as perturbative 
estimates are concerned these diagrams can be discarded, therefore, since they do not con- 
tribute to n\ growth of the coefficients. The remaining UV-divergent planar diagrams will 
contain among them diagrams containing chains of gluon self-energy insertions and other 
structures which must be combined with renormalon diagrams with chains of internal 
fermion bubbles and other structures to produce a gauge-invariant contribution propor- 
tional to a power of b, using the pinch technique or background field method |T2| . The 
planar diagrams of interest are those not involved in the construction of a gauge-invariant 
effective charge, therefore. The hope would be to understand why their contribution is 
'small'. This would not, unfortunately, explain why the RS-invariant effective charge 
beta-function coefficients are reproduced so well in the large- N limit, since this involves 
a combination of perturbative coefficients and beta-function coefficients. The piece that 
is so well approximated in the large- N limit is that which dominates in a Banks-Zaks 
expansion in e = — Nf) around e = [40| where asymptotic freedom is lost. This 
has been studied in the context of IR freezing of observables since for suitably small e 
there is expected to be an IR fixed point |37] . 
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5 Conclusions 



In this paper we have proposed an improvement of the renormalon-inspired 'leading-5' 
resummations of QCD perturbation theory which have been previously used by various 
authors [1-6] to assess the reliability of fixed-order perturbative predictions. Such re- 
summations are RS-dependent under the full QCD RG transformations. To avoid this 
difficulty the strategy is to approximate the RS-invariant effective charge beta-function 
coefficients by retaining their 'leading-6' part, which is completely determined by exact 
all-orders large- Nf results. Fixed-order perturbative approximations in any RS can then 
be obtained from the approximated RS-invariants by using the exact QCD RG. If the 
exact NNLO invariant is known it can be included. In this way the resummation includes 
the exact NLO and NNLO perturbative coefficients in any RS. 

The RS-invariant resummation was performed for the e + e~ i?-ratio, R T the analo- 
gous decay ratio for the tau-lepton, and DIS sum rules. Comparison with fixed-order 
perturbation theory in the effective charge RS revealed impressive convergence to the 
resummed result for the e + e~ i?-ratio at LEP/SLD energies, Q = 91 GeV. As the value 
of Q was reduced oscillatory behaviour of the fixed-order results above and below the 
resummed value was increasingly evident, reflecting the alternating-sign factorial growth 
of the perturbative coefficients resulting from the dominant UVi renormalon singularity. 
Even at Q — 1.5 GeV the resummed value was reasonably approximated until ninth order 
perturbation theory. 

For R T (Q = m T ), which is also UVi dominated, there was also a satisfactory approx- 
imation to the resummed value, although with much larger oscillations than for R at a 
comparable value of Q, and with an earlier breakdown of perturbation theory beyond 
fifth-order. 

In contrast DIS sum rules which have an IRi infra-red singularity exhibited much less 
satisfactory behaviour with successive orders moving steadily away from the resummed 
result, reflecting the fixed-sign factorial growth of the coefficients. 

Using the difference between the exact NNLO EC approximation and the resummation 
to estimate the uncertainty with which A-^g could be determined indicates that for R at 
Q = 91 GeV a s (M z ) could be determined to three significant figures with ideal data. 

For R T one concludes that d~a s (Mz) = 0.001 from the NNLO-resummed difference. 
This is a much smaller uncertainty than deduced by Neubert |J from a comparison with 
the naive RS-dependent leading-6 resummation. The RS-dependence means that the naive 
resummation is sensitive to the MS scheme a s {m T ) being assumed for the coupling. Other 
a priori reasonable choices would dramatically change the resummed result, and hence 
we would argue that this estimate of the uncertainty is too pessimistic. 

We regard the impressive performance of fixed-order QCD perturbation theory for the 
UV-renormalon dominated quantities as the key result of this analysis. 

Various technical issues related to the resummation and possibilities for future devel- 
opments were also discussed. 



21 



Acknowledgements 

We would like to thank David Barclay, Georges Grunberg, Jan Fischer, and Andrei Kataev 
for a number of stimulating discussions. 

D.G.T gratefully acknowledges receipt of a P.P.A.R.C. U.K. Studentship. 



22 



References 

C.N.Lovett-Turner and C.J.Maxwell, Nucl.Phys. B452 (1995) 188. 
M.Beneke and V.M.Braun, Phys.Lett. B348 (1995) 513. 
P.Bali, M.Beneke and V.M.Braun, Nucl.Phys. B452 (1995) 563. 
M.Neubert, CERN-TH 7524/94 ||hep-ph/9502"26j ; Phys.Rev. D51 (1995) 5924. 
M.Neubert, CERN-TH/95-112 ||hep-ph/950943l . 



M.Girone and M.Neubert, CERN-TH/95-275 [|hep-ph/9511392 
M.Beneke, Nucl.Phys. B405 (1993) 424. 
D.J.Broadhurst, Z.Phys. C58 (1993) 339. 
D.J.Broadhurst and A.L.Kataev, Phys.Lett. B315 (1993) 179. 
M.Beneke and V.M.Braun, Nucl. Phys. B426 (1994) 301. 
M.Neubert and C.Sachrajda, Nucl.Phys. B438 (1995) 235. 
N.J.Watson, Phys.Lett. B349 (1995) 155. 

C. N.Lovett-Turner and C.J.Maxwell, Nucl. Phys. B432 (1994) 147. 

K.G.Chetyrkin, A.L.Kataev and F.V.Tkachov, Phys.Lett. B85 (1979) 277; M.Dine 
and J.Sapirstein, Phys.Rev.Lett. 43 (1979) 668; W.Celmaster and R.J.Gonsalves, 
Phys.Lett. B44 (1980) 560. 

S.G.Gorishny, A.L.Kataev and S.A.Larin, Phys.Lett. B259 (1991) 144; 
L.R.Surguladze and M.A.Samuel, Phys.Rev.Lett. 66 (1991) 560; 66 (1991) 2416 (E). 

J.Chyla, Phys.Lett. B356 (1995) 341. 

S.G.Gorishny and S.A.Larin, Phys.Lett. B172 (1986) 109; E.B.Zijlstra and W. van 
Neerven, Phys.Lett. B297 (1992) 377. 

S.A.Larin and J.A.M.Vermaseren, Phys.Lett. B259 (1991) 345. 

D. T.Barclay, C.J.Maxwell and M.T.Reader, Phys.Rev D49 (1994) 3480. 
P.M.Stevenson, Phys.Rev. D23 (1981) 2916. 

G.Grunberg, Phys.Lett. B95 (1980) 70; Phys.Rev. D29 (1984) 2315. 
G.Grunberg, Phys.Rev. D46 (1992) 2228. 
M.Beneke, Phys.Lett. B307 (1993) 154. 
G.Grunberg, Phys.Lett. B325 (1994) 441. 



23 



LEP Collaborations Joint Report. CERN-PPE/93-157 (1993). 
R.Marshall, Z.Phys. C43 (1989) 595. 

G.Rahal-Callot, to appear in Proc.Int.Europhysics.Conf. on High Energy Physics, 
Brussels, Belgium, September 1995. 

S.R.Wasserbaech, to appear in Proc. Workshop on Tau/Charm Factory, Argonne, 
Illinois, June 1995. 

CCFR-NuTeV collaboration. Contribution to 30th Rencontres de Moriond: QCD 
and High Energy Hadronic Interactions, Meribel Les Allues, France, 19-25 March 
1995. ||hep-ex/95060l0 |. 



V.M.Braun and A.V.Kolesnichenko, Nucl.Phys. B283 (1987) 723. 

F. Le Diberder and A.Pich, Phys.Lett. B286 (1992) 147; F.Le Diberder and A.Pich, 
Phys.Lett. B289 (1992) 165. 

P.Raczka and A.Szymacha, Institute of Physics, Warsaw University, preprint 
IFT/1/95, ||hep-ph/9602245|1 . 

P.Raczka and A.Szymacha, Institute of Physics, Warsaw University, preprint 
IFT/13/94, ||hep-ph/9412"2"36] . 

G. Auberson and G.Mennessier, J.Math.Phys. 22 (1981) 2472; Commun.Math.Phys. 
100 (1985) 439. 

D.T.Barclay, private communication. 

A.C.Mattingly and P.M.Stevenson, Phys.Rev. D49 (1994) 437. 
P.M.Stevenson, Phys.Lett B331 (1994) 187. 

Yu.L.Dokshitzer, G.Marchesini and B.R.Webber. CERN-preprint CERN-TH/95-281 
frep-ph/ 95 12336 1. 



G.'t Hooft, Commun.Math.Phys. 86 (1982) 449. 
T.Banks and A.Zaks, Nucl. Phys. B206 (1982) 23. 



24 



